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Abstract. Anisotropic flow of hadrons is studied in heavy ion collisions at SPS and RHIC 
energies within the microscopic quark-gluon string model. The model was found to reproduce 

fT) \ correctly many of the flow features, e.g., the wiggle structure of direct flow of nucleons at midra- 
pidity, or centrality, rapidity, and transverse momentum dependences of elliptic flow. Further 
predictions are made. The differences in the development of the anisotropic flow components 

"P^ . are linked to the freeze-out conditions, which are quite different for baryons and mesons. 

Oh: 

p; 1. INTRODUCTION 

cu : 

The study of properties of extremely hot and dense nuclear matter, and the search for 
anticipated transition to a deconfmed phase of quarks and gluons, the so-called Quark- 
^ ■ Gluon Plasma (QGP), is one of the main objectives of heavy ion experiments at ultra- 
relativistic energies. Both theorists and experimentalists are looking for genuine QGP 
fingerprints, that cannot be masked or washed out by processes on a hadronic level. At 
present, the expansion of highly compressed nuclear mater in the direction perpendicular 
to the beam axis of the colliding heavy ions, known as collective flow, is believed to be 
one of the most promising signals to detect the creation of the QGP j 1 121.3 j . Since the 
development of flow is closely related to the equation of state (EOS) of nuclear matter, 
the investigation of the flow can shed light on the transition to the QGP phase accom- 
panied by its subsequent hadronization |4|5lb|7l8l9|lU|ll|12|13|14|15|lb117|18|19| . If the 
transition from the QGP to hadronic phase is of first order, the vanishing of the pressure 
gradients in the mixed phase leads to the so-called softening of the EOS |6|7j . The latter 
should be distinctly seen in the behavior of the excitation function of the collective flow. 
This circumstance explains the great interest in the transverse flow phenomenon. 

The Fourier expansion technique is usually employed to study collective flow phenomena 
since |2U|21| . The invariant distribution Ed 3 N/d 3 p is presented as 
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where pt and y are the transverse momentum and the rapidity, and <p is the azimuthal 
angle between the momentum of the particle and the reaction plane. The first two Fourier 
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coefficients in Eq. ((TJ), v\ and t>2, are dubbed directed flow and elliptic flow, respectively. 
Since both types of anisotropic flow depend on rapidity y, transverse momentum p t , and 
the impact parameter of an event b (i.e., v n = v n (xj), where {2^=1,2,3} = {y,Pt,b}), the 
following differential distributions are usually applied 

, a x f x f ] , , d 3 N , 2 / r*f> d 3 N , 2 
v n (xi,Ax,#) =J cos(n0) — -d Xfri / / -g-dx j+i . (2) 

Model calculations suggest that elliptic flow is built up at the early phase of nuclear 
collisions |10J14pi5j . whereas directed flow develops until the late stage of the reaction [221 
I23|24j . But it is well known that the particles with high transverse momentum are emitted 
at the onset of the collective expansion, i.e., their directed flow can carry information 
about the EOS of the dense nuclear phase. The study of the collective flow development 
is, therefore, closely connected to the freeze-out picture. In this article the microscopic 
quark-gluon string model (QGSM) |25|2fci| is employed to investigate the formation and 
evolution of anisotropic flow components in heavy ion collisions at SPS (^/s = 17.8 AGeV) 
and RHIC (\/s = 130 and 200 AGeV) energies. Note that QGSM does not implement 
the formation of a QGP at the early stage of the collision. Our goal is to understand to 
what extent the characteristic signals of the hot nuclear matter can be reproduced. A 
noticeable discrepancy between experimental data and QGSM predictions being observed, 
this should be considered as an indication for new processes not included in the model. 



2. MODEL 

The QGSM is based on the 1/N C (where N c is the number of quark colors or fla- 
vors) topological expansion of the amplitude for processes in quantum chromodynamics 
and string phenomenology of particle production in inelastic binary collisions of hadrons. 
The diagrams of various topology, which arose due to the 1/N C expansion, correspond 
at high energies to processes with exchange of Regge singularities in the t-channel. For 
instance, planar and cylindrical diagrams corresponds to the Reggeon and Pomeron ex- 
change, respectively. The QGSM treats the elementary hadronic interactions on the basis 
of the Gribov- Regge theory (GRT), similar to the dual parton model |2Zj and the VENUS 
model |2H]- This implies the consideration of subprocesses with quark annihilation and 
quark exchange, corresponding to Reggeon exchanges in two-particle amplitudes in the 
GRT, and with color exchange, corresponding to the one and more Pomeron exchanges 
in elastic amplitudes. The hh collision term includes also single and double diffraction 
subprocesses, antibaryon-baryon annihilation and elastic scattering, as well as the hard 
gluon-gluon scattering with large Q 2 > 1 (GeV/c) 2 momentum transfer [21]. 

The inelastic hh cross section cjj n (s) can be calculated via the real part of the eikonal 



2vr f {l - exp [-2u R (s, b)] } bdb . (3) 
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Here s is the center-of-mass energy of the reaction. The eikonal u(s, b) can be presented 
as a sum of three terms corresponding to soft and hard Pomeron exchange, and triple 
Pomeron exchange, which is responsible for the single diffraction process, 

u R (s, b) = u* oft (s, b) + «L(s, b) + uf nple (s, b) . (4) 
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Using the Abramovskii-Gribov-Kancheli (AGK) cutting rules jSUj the inelastic cross sec- 
tion of hh interaction can be presented as 

0~in{s) = a ijk(s) , (5) 

i,j,k=0;i+j+k>l 

oo 

o'ijk(s) = 2tc J bdbexp 



The last equation enables one to determine the number of cut soft and hard Pomerons, i.e., 
the number of strings and hard jets. The single Pomeron exchange leads to the formation 
of two quark-diquark or quark- ant iquark strings. With rising energy the processes with 
multi- Pomeron exchanges become more and more important. The contribution of the 
cylinder diagrams to the scattering amplitude increases like s ap ^°' > ~ 1 , while that of the so- 
called chain diagrams corresponding to n- Pomeron exchanges (n > 2) rises like s^H )- 1 ! 
with «p(0) > 1 being the intercept of a Pomeranchuk pole. Strings, which are formed 
in the course of a hh or A+A collision, decay later into secondary hadrons. Similar to 
hadronic collisions, string fragmentation into hadrons proceeds independently in A+A 
collision also. However, in the latter case these hadrons are allowed to interact with other 
hadrons after a certain formation time, while the valence quarks and diquarks can interact 
promptly with the reduced cross sections. The model simplifies the nuclear effects and 
concentrates on hadron rescattering. As independent degrees of freedom QGSM includes 
octet and nonet vector and pseudoscalar mesons, and octet and decuplet baryons, and 
their antiparticles. 

The transverse motion of hadrons in the QGSM arises from different sources: (i) pri- 
mordial transverse momentum of the constituent quarks, (ii) transverse momentum of 
(di)quark-anti(di) quark pairs acquired at string breakup, (iii) the transverse Fermi mo- 
tion of nucleons in colliding nuclei, and (iv) rescattering of secondaries. Parameters of the 
first two sources are fixed by comparison with hadronic data. The Fermi motion changes 
the effective transverse distribution of strings formed by the valence quarks and diquarks 
of the target and projectile nucleons. Thus, the original strings are not completely parallel 
to the beam axis. Further details of the QGSM can be found elsewhere |25|26j . We start 
from the study of energy and centrality dependence of the directed flow. 

3. DIRECTED FLOW 

The directed flow of nucleons at BEVALAC/SIS energies {E lah = 0.1 AGeV - 1 AGeV) 
and at AGS energy (Ei ab = 10.7 AGeV) has a characteristic S-shape attributed to the 
standard (p x /A) distribution. It grows linearly with rising rapidity between the target 
and projectile fragmentation regions. Conventionally, we will call this type of flow, for 
which the slope dv\/dy cm is positive, normal flow, in contrast to the antiflow for which 
dvi/dycm < in the midrapidity region. 

The one-fluid hydrodynamic models indicate that deviations from the straight line be- 
havior of the nucleon flow can be caused solely by the softening of the EOS due to the 
QGP creation in microscopic string model calculations such deviations were first ob- 
served in very peripheral Au+Au collisions at AGS energy |£L It appeared, however, that 
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the effect is shifted to more central topologies [23 24J as the collision energy increases. The 
phenomenon leading to the formation of a characteristic wiggle structure (3T] of the di- 
rected flow is caused by shadowing, which plays a decisive role in the competition between 
normal flow and antiflow in noncentral nuclear collisions at ultrarelativistic energies. 
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Figure 1. Directed flow of nucleons (full cir- 
cles) and pions (open circles) function of 
rapidity in lead-lead collisions at 160 AGeV. 
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Figure 2. Directed flow of protons and pi- 
ons in Pb+Pb collisions at 40 AGeV and 160 
AGeV measured by NA49 Collaboration 



Hadrons, emitted with small rapidities in the antiflow area, can propagate freely, while 
their counterparts will be absorbed by dense baryon-rich matter. The rapidity dependence 
of the directed flow of nucleons and pions in Pb+Pb at SPS is presented in Fig. ^for six 
different impact parameters. Here deviations from the straight line start to develop at 
b/b max = 0.45 in the central rapidity window. Recently, the wiggle structure of directed 
flow of protons in peripheral lead-lead collisions at SPS has been observed by NA49 
Collaboration [32] (see Fig. |2J). This result cannot be explained by the effective softening 
of the EOS in microscopic models due to production of strongly interacting string matter, 
because the number of the strings is significantly reduced as the collisions become more 
peripheral. In hydrodynamics, even if the tilted region is non-perpendicular to the beam 
axis ^7] , the tendency should be opposite: Here the antiflow reaches its maximum strength 
in semicentral events and almost disappear in peripheral ones. 

Our next step is to study the time development of the directed flow of strange hadrons, 
primarily, kaons and antikaons, in order to check the role of baryon-antibaryon asymmetry. 
For instance, K~ and K can be absorbed via the channels such as K~ +p — > A + 7r + , etc., 
whereas there are no analogous reactions for K + and K°. How important is this reaction 
asymmetry at SPS and, especially, at RHIC energies, where the matter is expected to be 
meson-dominated? These problems have been addressed in [33J. 

The time evolution of directed flow of kaons and antikaons in minimum bias Pb+Pb 
collisions at 160 AGeV is presented in Fig. El Here the coefficient vf{y) is calculated 
in different transverse momentum intervals at early, t = 3 fm/c and t = 10 fm/c, and 
the final stage of the reaction, t > 60 fm/c. To avoid ambiguities, all resonances in the 
scenario with early freeze-out were allowed to decay according to their branching ratios. 
The flow evolution is seen quite distinctly. At early stages of the collision directed flow 
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of both kaons and antikaons is oriented in the direction of normal flow similar to that 
of nucleons (231 • Within the error bars there is no differences between (K + + K°) and 
(K~ + K ). At this stage the matter is quite dense, mean free paths of particles are short, 
and similarities in kaon production and rescattering dominate over inequalities caused by 
different interaction cross-sections. It is worth mentioning that the directed flow of K and 
K is already sizable at t = 3 fm/c. This can be explained by a kick-off effect associated 
with the early stage of the collision, when the nuclei pass through each other. Later on the 
system becomes more dilute. For both kaons and antikaons the directed flow experiences 
significant transformations. Already at t — 10 fm/c the antiflow of antikaons starts to 




Figure 3. (y) (solid circles) and t>f (y) 
(open circles) in min. bias Pb+Pb collisions 
at SPS energy in high p t interval 0.6 < pt < 
0.9 GeV/c (left) and for all pt (right panels) 
at times t = 3 fm/c, 10 fm/c, and final. 
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Figure 4. The same as Fig. 3 but for mini- 
mum bias Au+Au collisions at RHIC (y/s = 
130 AGeV). 



built up in the midrapidity range. Note, that t = 10 fm/c corresponds to the maximum 
of the kaon dN/dt distribution over their last elastic or inelastic interaction [33]. Here the 
differences in interaction cross sections and possible reaction mechanisms become crucial. 
Due to larger interaction cross-sections of K~ and K with other hadrons, the directed 
flow of these particles changes the orientation from a weak normal to strong antiflow. 
Even (K~ + K ) with high transverse momentum demonstrate distinct antiflow, while 
the flow of (K + + K°) remains almost unchanged compared to that at t — 10 fm/c. 

The directed flow of kaons and antikaons in minimum bias Au+Au collisions at \fs = 
130 AGeV is displayed in Fig. |U again at early stages, t = 3 fm/c and t = 10 fm/c, 
and at the final one, t > 100 fm/c. It is interesting that at t — 3 fm/c (i) the flow of 
(K + + K°) coincides within the statistical errors with the (K~ + K°) flow, and (ii) the 
flow is generally very similar to that at the SPS energy at time t = 3 fm/c. Except of the 
target and projectile fragmentation region, where again the flow is probably produced by 
the initial kick, the kaon flow at this early stage of gold-gold collisions at RHIC energy 
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is isotropic with respect to the impact parameter axis. The spatial anisotropy in the 
distribution of baryonic charge seems to be unimportant at this stage. At t = 10 fm/c 
not only the directed flow of antikaons, but also that of kaons becomes antiflow-aligned at 
midrapidity. Similar behavior has been found within the RQMD model for the directed 
flow of nucleons at RHIC [ST] , suggesting that the nucleon directed flow is a side effect of 
the elliptic flow. The flow of produced particles, pions |31|35| and kaons [S3], was found 
to be very flat at \y\ < 2, in stark contrast to the QGSM predictions. We are awaiting 
the experimental data to resolve this problem. 



4. ELLIPTIC FLOW 



Rapidity and pseudorapidity distributions of elliptic flow of pions and charged particles 
at RHIC |I5J are depicted in Fig. For both energies elliptic flow displays strong in-plane 
alignment in accordance with the predictions of Ref. [S]. In the mid(pseudo)rapidity the 
flow is almost constant. It rises slightly at \y\, \r]\ « 1.7, and then drops with increasing 
rapidity. Pseudorapidity dependences of the elliptic flow of charged particles in the whole rj 
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Figure 5. Rapidity dependence of elliptic 
flow of pions and charged particles at RHIC 
energies. See text for details. 



Figure 6. Centrality dependence of elliptic 
flow of pions and charged particles at RHIC 
energies. 



range, which were obtained before the experimental data ^H] , are in a good agreement with 
the the results reported by PHOBOS Collaboration 36 . To elaborate on the influence 
of hard processes and multi-Pomeron exchanges on the elliptic flow formation the flow 
caused by the subprocesses without the hard and multichain contributions is also plotted 
in Fig. E| It seems that in Au+Au collisions at yfs = 130A GeV the magnitude of 
the signal (except of the midrapidity range) can be reproduced without the many-string 
processes, although the particle multiplicity in the latter case is reduced by 30%. At 
■\/s = 200A GeV their role becomes more significant, because the elliptic flow caused by 
other subprocesses cannot exceed the limit of 3—3.5% . Here it is important to stress 
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that the multichain diagrams alone, without rescattering, cannot affect the elliptic flow 
at all. The flow increases solely due to secondary interactions (absorption or rescattering) 
of produced particles in spatially asymmetric systems. 

Figure El presents the centrality dependence of elliptic flow of charged particles. Since 
the centrality of events in the experiment [37| has been determined via the ratio of charged 
particle multiplicity to its maximum value N c h/ (N c h) max , we compare the Vi \N c hj {N c h) max ] 
signal with the original impact parameter dependence V2Q)). One can see that as a func- 
tion of the impact parameter b elliptic flow is saturated at b « 8 fm for both energies, 
while as a function of the multiplicity ratio it increases nearly linearly with decreasing 
multiplicity up to N c h/(N c h) max ~ 0.15. The agreement with the STAR Collaboration 
data, obtained by 2-cumulant and 4-cumulant method |37|38| . is good; similar results for 
the centrality dependence of t>2 were obtained also by the PHENIX Collaboration |3H|. As 
expected, the flow in the midrapidity region is caused mainly by pions. The magnitude of 
the pionic flow in the QGSM calculations is twice as large as obtained, e.g. with RQMD 
|4Uj . Without many-string processes the QGSM is able to describe the flow only in central 
and semicentral collisions. It predicts a drop of the elliptic flow as the reaction becomes 
more peripheral, which is similar to the predictions of other string models. Note that 
neither the pseudorapidity nor the centrality dependence of the elliptic flow at RHIC is 
reproduced correctly so far by the hydrodynamic models. 

Figure [7| depicts the pr- dependence of the 
elliptic flow of strange hadrons at RHIC. 
The flow is close to zero for hadrons with 
Pt < 0.25 GeV/c, then rises linearly 
(kaons) up to v^^Pt) ~ 10% within the 
interval 0.25 < px < 1.5 GeV/c, and sat- 
urates at p T > 1.5GeV/c in accord with 
the experimental data jH]. The flow of 
lambdas is weaker than the kaon flow at 
Pt < 1-3 GeV/c. At higher transverse mo- 
menta the situation is changed: here the A 
flow is stronger. This is the general trend 
in the p? dependences of mesonic and 
baryonic flow at RHIC energies. The rel- 
ative reduction of the elliptic flow within 
the interval 1 < px < 2 GeV/c can be ex- 
plained by the interplay between the flow 
of high-p r particles, emitted at the onset 

of the collision, and the hydro-type flow of particles, which gained their transverse mo- 
mentum in subsequent secondary interactions. Therefore, for proper understanding of the 
flow evolution it is very important to investigate the freeze-out picture of the hadrons. 

5. FLOW AND FREEZE-OUT 

The dN/dt distributions of n c h, tt, N, A, which are decoupled from the fireball after the 
last elastic or inelastic collision, are shown in Fig. |H1 For the analysis ca. 20000 gold-gold 
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collisions with the impact parameter b = 8 fm at y/s = 130 AGeV were simulated. One 
can see, that the particles are emitted during the whole course of the system evolution. 
In this respect the freeze-out picture obtained in microscopic model is different from the 
sharp freeze-out assumed in hydrodynamic models. Compared to AGS and SPS energies 
[3*4"] . a substantial part of hadrons leave the system immediately after their production 
within the first two fm/c. The pion distribution has a peak at t = 5 — 6 fm/c, while the 
distributions of baryons are wider due to the large number of rescatterings, which shift 
their dN/dt maxima to later times, t = 10 — 12 fm/c. Since the pion fraction dominates 
the other hadrons, the distribution of charged particles is similar to the pion one. 
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Figure 8. dN/dt distribution of n c h,Tr,N,A 
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and elliptic flow of these particles (right) for 
Au+Au collisions with b = 8 fm at yfs = 130 
AGeV. 



Figure 9. Gold-gold collisions with 6 = 8 
fm at y/s = 130 AGeV: (a) Snapshots of the 
v~2 at 2, 4, ... 80 fm/c (from bottom to top), 
and (b) elliptic flow of pions already frozen 
at 2, 4, ... 20 fm/c (from top to bottom). 



The elliptic flow carried by these species is presented in Fig. |H1 as well. The baryonic 
and mesonic components are completely different: pions emitted from the surface of the 
expanding fireball within the first few fm/c carry the strongest flow, while later on the flow 
of pions is significantly reduced. In contrast to this, the baryon fraction acquires stronger 
elliptic flow during the subsequent rescatterings, thus developing the hydro-like flow. The 
saturation of the flow at the late stages can be explained by the lack of rescattering, since 
the expanding system becomes more dilute. 

Finally, the time evolution of the elliptic flow of pions at midrapidity is studied. Two 
varieties of rapidity distribution are displayed in Fig. El elliptic flow of pions frozen at t = 2 
fm/c, 4 fm/c, etc. till the late stages of the reaction; and the rapidity profile of the total 
pion elliptic flow at the same times. In the latter case, similarly to Sec. El the interactions 
were switched off, but resonances were allowed to decay. Here two observations should be 
mentioned. Firstly, although the elliptic flow at y = reaches the maximum already at t ~ 
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6 fm/c, its formation is not over due to continuous freeze-out of particles. Secondly, there 
is no one-to-one correspondence between the apparent elliptic flow and the contribution 
to the final flow coming from the "survived" fraction of particles. For instance, apparent 
elliptic flow of pions at t — 2 fm/c is weak, but pions which are frozen at this moment 
have the strongest elliptic anisotropy caused by the absorption of the pion component in 
the squeeze-out direction. 

6. CONCLUSIONS 

In summary, the features of directed flow in the model can be stated as follows: 

(i) wiggle structure of nucleon flow at midrapidity (SPS, RHIC); 

(ii) the negative slope of the flow at \y\ < 2 similar for all hadrons (RHIC); 
(hi) normal flow of high-p t hadrons at midrapidity (SPS, RHIC); 

(iv) flows of kaons and antikaons are different at SPS but similar at RHIC. 
Elliptic flow features at RHIC: 

(i) flow &ty/s = 200 AGeV is only 10-15% stronger than that at y/s = 130 AGeV; 

(ii) V2(rj) varies slightly within the interval \r]\ < 2 and quickly drops at \r]\ > 2; 

(hi) V2(pt) of mesons is stronger than the baryonic flow at p t < 1.5 GeV/c but saturates 
earlier; 

(iv) the time evolution of the mesonic flow is opposite to that of the baryonic one. 
Many of these peculiarities are linked to different freeze-out pictures for baryons and 
mesons. In particular, it is shown that the flows of particles emitted at the beginning of the 
collision and during the course of the reaction are quite different. Therefore, phenomena 
such as collective flow and particle freeze-out should not be considered independently. 
The collective flow of hadrons appears to have a multi-component structure, caused by 
rescattering and absorption in a spatially anisotropic medium, which deserves further 
investigations. 
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